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A new type of approximate scaling compatible with the Balitsky-Kovchegov equation with run- 
ning coupling is found, which is different from the previously known running coupling geometric 
scaling. The corresponding asymptotic traveling wave solution is derived. Although featuring dif- 
ferent scaling behaviors, the two solutions are complementary approximations of the same universal 
solution, and they become equivalent in the high energy limit. The new type of scaling is observed 
in the small- a; DIS data. 



I. INTRODUCTION 

Due to intense soft gluon radiation, hadronic processes at high energy involves strong color fields which have a 
nonlinear evolution, with in particular saturation effects. The mixing between quantum and nonlinear effects makes 
^v^j this high energy behavior difficult to derive from the QCD Lagrangian. However, in the last decade, many progresses 
have been made, in particular for the case of dense-dilute collisions, such as deep inelastic scattering (DIS) or proton- 
nucleus collisions. Two dual formulations of the same renormalization group evolution with rapidity for such collisions 
have been built, i.e. Balitsky's hierarchy [l[ and the JIMWLK equation @, H, 0, H, @> 0] • These equations resum, 
at the leading logarithmic (LL) accuracy, the soft gluon emissions, taking into account nonlinear high-density effects 
at all orders for one of the two colliding particles. Neglecting correlations in the dense regime, those two formalisms 
reduce to the Balitsky-Kovchegov (BK) equation [H, H, [9[ . In fact, numerical simulations of the JIMWLK equation 
' show no significant differences from solutions of the BK equation [l(|. The BK equation gives the evolution with 
rapidity of the dipole-target amplitude, appearing in the dipole factorization of DIS, for example. The geometric 
Q_T scaling [ll| property, 

& <Jdis{Q\ Y) ~ a (s f (\og(Q 2 /A 2 QCD ),Y)) (1) 

^! with s f (L,Y) = L — avY , (2) 

I \ seen in the DIS data at high rapidity, can be explained by the scaling N(L, Y) = N(sf(L, Y)) of the weighted Fourier 
transform Q of the dipole-target amplitude, which is approximately verified by the generic solutions of the BK and 
' of the Balitsky- JIMWLK equations [12 ], We use the notation a = N c a s /ir, and L = log(kt 2 /AQ C£) ), kt being the 
, Fourier conjugate of the parent dipole size r. 

In the three last years, a lot of efforts have been devoted to the improvement of the saturation equations, in m ainly 
two directions. First, the calculation of next-to-leading logarithmic (NLL) contributions has started [ll, [HI, lla j. 
which are certainly required to achieve precision on the theory side. Moreover, at the LL accuracy, the QCD coupling 
remains fixed. Hence, the BK equation at LL accuracy is not expected to give good results over a wide energy 
range, due to the lack of asymptotic freedom. NLL contributions provides however indications about the best way to 
implement running coupling in the saturation equations [HI, [TiJ . Running coupling effects significantly modifies the 
whole behavior of the solutions [12, HH Iffi [T^ , whereas NLL contributions to the kernel or to the nonlinear term 
become irrelevant at high enough energy [19(. 

Second, the study of the inclusion of Pomeron loops in the saturation equations, and of the fluctuation effects 
induced by them, has been very active. However, it has been recently realized that in the case of dense-dilute 
processes, running coupling effects reduce the fluctuation effects in the high energy limit [20ll2l|. and even kill them 
completely at phenomenologically relevant rapidities [20|. Therefore we won't discuss further those effects in the 
present paper. 

For the case of inclusive DIS, one has to calculate the forward part of the dipole-target amplitude. Assuming impact 
parameter independance or a factorized momentum-transfert dependence [22], the BK equation for the forward 
dipole-target amplitude writes 

d Y N(L,Y) = a X (-d L )N(L,Y)-aN 2 . (3) 
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In the previous formula, x(t) = 2\P(1)— ^(7)— \&(1— 7) is the characteristic function of the BFKL kernel. Considering 
the paren 
coupling 



the parent dipole relative kt as the relevant scale for the coupling, one obtains the forward BK equation with running 



d Y N(L, Y) = ±. X (-8 L )N(L, Y) - ^N 2 , (4) 

where b = (11N C - 2N f )/12N c . 

Our main tool to study Eqs. (|3I4|) is the traveling wave method [H, [24|, borrowed from nonlinear physics. The BK 
equation ([3J shares indeed the same qualitative structure with the FKPP equation (2f| [26| , for which that method 
has been developped, with an instability of the vacuum, a diffusion term, and a nonlinear damping. Thanks to this, 
those two equations have both traveling wave solutions fl8l H3] , with the following property. For any initial condition 
with a steep enough tail (which is the case in QCD due to color transparency) the traveling wave solution will evolve 
to an universal asymptotic solution, loosing memory of the initial condition. This universality property constrains 
not only the limit of the solution, but also the behavior of the solution at the first two subleading orders. This whole 
universal behavior is known [28| for the BK equation with fixed coupling ([3]). In the running coupling case, similar 
universal asymptotic traveling wave solutions appears, although the convergence towards the asymptotic behavior is 
less spectacular. The solution has the approximate scaling property N(L, Y) = N(s g (L, Y)) [l2j with the variable 

fuY 

s g (L,Y) = L-^j—, (5) 

for large L and s g <C VY, which becomes exact in the large L and large Y limit. This corresponds to a geometric- 
scaling property with log Q 2 ~ VY, instead of log Q 2 S oc Y for the fixed coupling geometric- scaling Q . In the running 
coupling case, the universal asymptotic solution is known including the first subleading order [161 1 1 81] only 1 . 

From the theory point of view, the BK equation with running coupling is the most efficient tool to derive analytically 
the dipole-target amplitude. However, numerical simulations [23 . |30j of the BK equation with running coupling have 
pointed out that the scaling function is significantly steeper near the saturation scale in the running coupling case 
than in the fixed coupling case, whereas the running coupling solution from [l(| [l8[ and the fixed coupling solution 
have similar shapes, with in particular the same coefficient j c ~ 0.62 playing the role of an anomalous dimension. 
That fact seem to rise some doubts about the relevance of the previously known asymptotic solution of ((4|): either 
this solution is not the same one as in the simulations, either strong corrections have been missed. Moreover, the 
running coupling geometric scaling (O has been shown compatible with DIS data [3lj . but not significantly better 
than the fixed coupling geometric scaling ^ . The too narrow Q 2 range of the available low- x data can explain partly 
that observation, but maybe not completely, as fixed coupling and running coupling solutions differ significantly. The 
purpose of the present paper is then to reexamine the problem of finding analytical approximate solutions of the BK 
equation with running coupling and understanding their validity range. 

For sake of simplicity, let us give the main results of our study, whose complete derivation is rather technical and 
is given in the text. We perform a search for scaling variables compatible with the equation Q. As a main result we 
find that, contrary to the fixed coupling case ([3]), no exact scaling can be obtained, but only approximative ones. We 
find two different versions of asymptotic scaling valid in the vicinity of the saturation scale and above, including the 
previously known one ([5]) and a new one 

for large L and s n <C L, which is not, rigorously speaking a geometric- scaling variable since it cannot be rewritten 
(by exponentiation) as kt 2 /Q 2 (Y). In Eqs. (|2|5l6p . the index / stands for fixed coupling, the index g for geometric 
and the index n for new. The as ymp totic traveling wave solution associated with the scaling variable ((6|) is derived 
using the traveling- wave method [23l . |24| . It writes 



N(L,Y) oc e -7c^+o(iogL) 



Ai Ui + 



(DL)V3 



(7) 



with -s n = l- V ^-%D L) V*, (8) 



However, some unpublished preliminary results concerning the second subleading order have been obtained by Stephane Munier. 



for large L and l<s„< \J~L. The neglected terms . . . are of order 0(L" 1 l' i ) when s n — O (L 1 / 3 ). The parameters 
7c v c , and D are defined in the text (see Eqs. (|44l45|A8p ). and Ai is the Airy function. The two asymptotic traveling 
wave solutions, associated with the scaling variables s g and s n are two complementary approximations of the same 
exact solution. In particular, they both include the universal features of the exact solution, but differ significantly 
concerning some non universal features, such as the shape of the wave front near the saturation scale. Having two 
different approximate solutions allows a better control on the validity range of each, and can help to estimate a part 
of the errors. 

The plan of our analysis is the following. In section[llJ a search for scaling variables compatible with the BK equation 
with running coupling ([U is performed. In section llTTI the asymptotic traveling wave solution (|7l8p associated with the 
scaling variable © is derived using the traveling wave method [23l . [24j | . The approximate asymptotic traveling wave 
solutions associated with the scaling variables s g ((5|) and s n ([6]) are compared in section Hvl The phenomenological 
relevance of the new type of scaling s n is shown in section [Vj The last section contains the conclusion and an outlook. 



II. SCALING VARIABLES FROM SATURATION WITH RUNNING COUPLING 
A. Preliminaries: Scaling for a medium-dependent FKPP equation 

The FKPP equation has been useful [13] to understand the scaling properties of the BK equation at fixed coupling. 
Therefore, before considering the BK equation with running coupling Q, let us study the following modified FKPP 
equation, mimicking running coupling effects 

d t F{x, t) = - [d 2 x F{x, t) + F(x, t) - F 2 {x, t)] , (9) 
x 

where the operator d 2 +I plays the same role as the kernel %(— &l) in (J4|. Note that this equation might be interesting 
also as a mean field description of a reaction diffusion process in statistical physics when a particular medium implies 
a time dependence decreasing with the distance x. 

F is a scaling solution associated with the scaling variable s(x, t) if it writes 

F(x,t) =F s (s{x,t)). (10) 

Inserting the condition (fTQ]) in ([9|), one gets 

x(d t s(x,t))F^s) = (d 2 x s(x,t))F^s) + (d xS (x,t)) 2 F>>(s) + F s (s) - F?(s) . (11) 

As F s (s) — F 2 (s), Fi.(s) and F""(s) are scaling functions, i.e. depend only on s, the coefficients of these functions 
must be scaling functions, in such a way that F(x,t) = F s (s) remains really an exact scaling solution of @. All in 
all, finding an exact scaling variable s associated with an exact scaling solution (fTUf of © is equivalent to solving the 
equations 

xd t s(x, t) = h{s) (12) 
d x s(x,t) - Ms), (13) 

/i and f% being some arbitrary functions of s. 

In order to have a smooth scaling variable s(x, t), one has also to impose that 

d t d x s(x, t) = d x d t s(x, t) . (14) 

Eqs. (|12I13P give 

d t d x s(x,t) = - h(s) &(s) (15) 

X 

d x d t s(x,t) = - fi(s) f 2 (s) - \ h{s) . (16) 
X x z 

One concludes that the inhomogeneous FKPP equation © has no nontrivial exact scaling solution. 

Note that for the standard FKPP equation, the same approach would give the scaling conditions (|12ll3p without 
the x prefactor in (fl"2"|) . Hence, the two conditions would be compatible, and give the analog of the geometric scaling 
variable with \ogQ 2 (Y) cx Y and nothing else. 



Large x approximate scaling solutions are interesting. For x large, one can drop the last term in (| 16[) without 
assuming /i(s) = 0. Matching the 1/x terms in (fT5|) and (|16|) gives 

m _ m 

his) h{s) ' [ n 

then fi(s) and his) are proportional, which allows us to write 

v 

x dts(x,t) = — — hi s ) with v any real number (18) 

d x s[x,t) = his)- (19) 

h reflects the invariance of the Ansatz (I10p under reparamctrizations. One can choose a specific h function without 
loss of generality. Choosing hi s ) = 1 5 one has 

v 

x dts(x,t) = — — with v any real number (20) 
d x s{x,t) = 1. (21) 

As discussed before, the equation §§§ has no exact scaling solutions. Then, the equations ((2"0"|) and (|2ip cannot be 
solved simultaneously beyond a large cc approximation. 

Consider the following alternative 

i) An exact solution for (|21[) and an approximative one for (|20[) : One gets from (|21[) s(x,t) = x — <^(t), then 

x a tS (x, i) = -x <t>'it) = -4>'it) 4>it) (l + fM) . (22) 
Hence, the equation (|20|) is approximately verified if 

That equation can be solved for s <C </>(i), as 



= Vu(t-to), (24) 
io being an integration constant. Thus, we find the scaling 



s{x,t)=x-jv(t-to), (25) 

which is valid for x large and six, t) -C x. The second condition means that t — to = 0(x 2 ), hence this scaling can 
be valid only if x and t are both large. This solution corresponds, in the case of the BK equation with running 
cou pling, t o the solution with a modified geometric scaling involving the square root of the rapidity, as proposed in 

ii) An exact solution for (|20| and an approximative one for f|2 1|) : One gets from (|20p 



2.r 

Then, the equation ([21]) rewrites 



vt 

s(x,t) = (p(x) - — . (26) 



// s vt . ip(x) s(x,t) 

1 ~ tp'ix) + — 2 = tp'ix) + ^ - ^f 1 ■ (27) 

As previously, one finds an approximate solution for x large and s(x,i) -C x, which is 

to being an integration constant. Hence, for x large and s(x, ()<! we have the approximate scaling variable 

-(*.*) = 1-^^- (29) 



B. QCD case: the Balitsky-Kovchegov equation with running coupling 



Let us now consider the BK equation with running coupling Q with the same approach. In order to search for 
scaling solutions 

N{L,Y) = N.{a{L,Y)) , (30) 

it is convenient to write the kernel of (J1J) as 



oo 1 

X(-^)=E-X W W with 7 G]0,1[. (31) 



p=0 1 

Thus, if (— 7 — d£)s{L, Y) is a scaling function, then (— 7 — dL) p N s (s(L, Y)) is a scaling function, and then 
x(~ ®l)N s (s(L, Y)) is a scaling function too. Hence, the constraints 

bLdYs{L,Y) = h{s) (32) 
d L s(L,Y) = f 2 (s) (33) 

are sufficient conditions to have a scaling solution N s (s(L, Y)) of the equation ([4}, and they are likely to be necessary 
too. As the conditions (|3"2"j) and (|3"3")) are similar to (TT2"|) and (fT3j) . the results found for the equation © apply for the 
Balitsky-Kovchegov equation ([!]). 

We can thus give the following conclusions. 

• On the contrary to the fixed coupling case, the BK equation with running coupling seems to have no exact 
scaling solutions. But in fact it is not really a problem, as the travelling wave method requires only a family 
of asymptotic scaling solutions. Moreover, an exact scaling would depend on the behavior of the coupling at 
transverse scales of the order of Aqcd or smaller. As the BK equation has been derived in perturbation theory, 
it should not give such non perturbative information. Thus, even if exact scaling solutions of ((H) exist, beyond 
the large L approximation, they would be unphysical. 

• The BK equation with running coupling ^ has solutions with large L scaling similar to to the solutions (|25l29p 
of the modified FKPP equation |9]). The region near the front, where Y — 0(L 2 ), is compatible with two types 
of scaling behavior at large L, either the geometric scaling depending on the square root of the rapidity 



as derived in [12j, |lfj, ll7|, Ilia] , or a new non-geometric scaling 

L v(Y - Yp) 
Sn ~ 2 2bL 



(35) 



In both cases, the parameter Yq reflects subasymptotic contributions. We will sometimes drop it in the following of 
the present theoretical study, but it is relevant for phenomenological studies of scaling [3l|, [HJ . 



III. TRAVELING WAVE SOLUTION WITH THE NEW SCALING VARIABLE 

In this section we apply the traveling wave method [23|, H3] to derive the asymptotic solution of the BK equation 
with running coupling based on the scaling (|35|) at large L and Y = 0(L 2 ). The most technical part of the calculation 
is given in the Appendix [X] The calculation is similar to the ones done in (l8| , in particular for the running coupling 
solution based on the scaling (jMf . However, we use a slightly different presentation of the method, doing the 
calculation step by step, instead of taking an Ansatz. It is useful in order to clarify the validity range of the solution, 
and to understand the dynamics hidden in the Ansatz. 

The first step is to rewrite the equation (£[]) in term of the scaling variable s n . In order to find the solution with 
both the scaling and the subleading scaling violations, one has to keep a second variable. It is possible to keep the 
effective time t, following [IB, G3, EB| j except that now t = a{L)(Y — Yq) = (Y— Yo)/bL. However, it seems more 
interesting to keep L. On the one hand, one can verify that, in that case, this choice will drive us to neglect less 



terms than the choice of t, giving hope to obtain a more accurate solution. On the other hand, it allows us to have 
an easier control on the validity range of our result, as our hrst approximation is the validity of perturbation theory 
i.e. bL ^> 1. Note that for all positive L and Y > Yq, s n (L,Y) < L/2, and that s n — L/2 corresponds to Y = Yq. 
Finally, the equation (f?)) rewrites in terms of s n and L 

- ~ d Sn N(s n , L) = X (- (l- y) 9 *« - dL ) ^( s ™' L ) " ^ 2 ( s «' L ) ( 36 ) 

where N(s n (L, Y),L) = N(L, Y). 

In the traveling wave formalism for nonlinear equations, one distinguishes different phase-space regions depending 
on the dynamical mechanism at work. Let us examine two important such regions where it is possible to make 
analytical predictions on the form of the solution. These two regions are the "front interior" and the "leading edge" , 
following the standard notations of Ref. [2J] . 



A. "Front interior" of the traveling wave solution: 

As discussed in section UH the scaling in s n is valid only at large L. Let us first study the "front interior" of the 
solution, i.e. the region defined by s n = 0(1). The solution of (|36p here admits a large L expansion in the front 
interior 



Inserting (pTFf in equation (|36p. one finds that A/q(s„) verifies 

- | ds n N {s n ) = X (-d Sn )N (s n ) - N^s n ) . (38) 

This equation can be solved numerically for any relevant value of v. Nevertheless we will need later an analytic 
expression for the tail of No(s n ) in the dilute regime, i.e. the solution of 

- \ d Sn N (s n )~x(-d Sn )N (s n ). (39) 
Factorizing the leading exponential decay (with some unknown parameter 7), one writes 

No(s„) = e"T s » f FI (s n ) , (40) 
fFi(s n ) behaving as a power for s„ large. The equation (|39|) then becomes 

^ fFl(Sn) ~ ~ d s f FI ( Sn )= X (l) fFl{s n ) - X '(7) d s jFl{s n )+Y,X {P \lV—r d P s jFl(Sn). (41) 

2 2 ^ p\ 



v= 2 -^l. (42) 

7 



One gets from (|4Tj) the dispersion relation 



Then, either 7 is arbitrary and all the derivatives of fpi have to vanish, either the first derivative is the only non- 
vanishing one, which gives 

v = 2 x'(l) (43) 

and thus 7 = j c , defined by 

Xilc) - 7c X'ilc) ■ (44) 

The velocity of the wave corresponding to "f c is 

„ = ^S). (45) 



Hence, the solutions of the equation (j3"9")l are 



N (s n ) = 4e^ s » with 7 ^ 7c (46) 
or N (s n ) = A (s„ + -B) e~ 7cS " , (47) 

with A positive and -B real. The solutions (|46|) are valid large L solutions of (f5B|) for s n <C \f~L. When s n — 0{\/~L) 
the terms with f' FI come at the same order as some term suppressed with s n /L, and the dynamics is modified. The 
solutions (f4"T|) are the large L critical solutions of ([3"6"|) . They are valid for s n <C L 1 / 3 . New dynamical effects occur 
for 0(s n ) = OIL 1 / 3 ), when the term with f' FI comes at the same order as a term suppressed with s n /L. This region 
is called the "leading edge" [11] and will be studied in the next subsection. 

Eqs. (|46l47p give both valid solutions of the linear part of the BK equation with running coupling However, if 
the initial condition is steep enough, the nonlinear damping term will select dynamically [23 , |24| ] the critical solution 
(jTHl , and then the contributions of the non-critical solutions ([4"6"|) will not survive asymptotically. In that case, the 
resulting solution is called a pulled front. By contrast, pushed fronts are solutions in which the generic partial waves 
of the type (|46p would survive, for exemple if the initial condition were not steep or if the phase velocity given by the 
dispersion relation would have no minimum. 

Going back to the equations (I36|) and (|37|) . one can in principle calculate Ni(s n ) and higher order terms recursively. 
But we will stop here in this work for two reasons. First, Ni(s n ) is not needed for the following calculations at the 
accuracy considered in this paper. Second, No{s n ) is free from NLL corrections (beyond the resummed corrections 
giving the running coupling considered here), whereas Ni(s n ) would depend on NLL contributions, N2(s n ) on NNLL 
contributions, etc. 



B. "Leading edge" of the traveling wave solution: 



We study now the equation (|36p when L and s n are both large. We drop the nonlinear term, because s n large 
corresponds to the tail of the front. Let us factorize the exponential decay of the solution as in (|40j) . but for the whole 
solution 



N(s n ,L) = e-^ s " f LE (s n ,L). 
Then, expanding the kernel around 7, the equation (|36p becomes 



(48) 



^ fLE(Sn,L) - \ 8 s J LE (s n ,L) = ]T ^X iP) h) 



p=0 



-(1-^)8.. -a 



7£n 

L 



/l_e(s„,L) 



(49) 



Hence for s n L, the cancelation of the leading term on each side gives again the dispersion relation (|42|) . and it 
remains 



P =i 



pi 



-11- 



8« - 8 L - 



/lb(s„,L) . 



(50) 



As we are interested by the solution selected by the front formation mechanism, which corresponds to (|47p . 9 s „/l_e 
has to be non-zero. Then, for f^E x s n /L 9 s „/lb, which corresponds a priori to s n <C VL, the leading terms in 
the equation ([50| gives 



2 8 Si J LE {sn,L) = -x'{l) 9 s „/lb(s„,L) 



(51) 



Hence, the relation (|4*3"|) holds while s„ \f~L, which fixes 7 = 7 C and v — v c = 2x( 7c )/7c- Then, the equation ([5 
becomes 



= x'(7c) 



~L - dL -~ 



1 



f LE ( Sn ,L)+J2-X {p) (lc) 



p=2 



pi 



L 



d s „ - d L 



IcSr 



For s n \[L, the leading terms in (|5^| gives 



X'(7c) 



7c Sr, 



L 



fL E (s n ,L) = - x"(7o) d 2 s J LE (s ni L). 



f LE (s n ,L). (52) 



(53) 



Hence, non-trivial effects appears for s n = ©(L 1 / 3 ), and Jle is allowed to have an i-dependance constrained by 



d L Jle < 0(L- 2 / 3 f LE ) 



(54) 



Thus, the leading term in f^E can be factorized as a prefactor exp(— times a scaling function depending 
only on Sn/L 1 / 3 , where (3 is a constant which have to be determined. This prefactor amounts to give a subleading 
correction to the traveling-wave scaling. We are thus led to define a more accurate scaling variable 



3/3L = 



L 
2 



v c Y 
2bL 



(55) 



The rather technical derivation of the leading edge solution is given in Appendix 1. 

Finally, the universal asymptotic traveling wave solution of the BK equation with running coupling associated with 
the scaling (l35|) writes, for large L and 1 <C s n <C \[L 



N(L, Y) cx e -icS„+o(io g L) Ai / ^ 



(DL) 1 / 3 



with 



L_v_cY_3b 1/3 
2 2bL 4 v ' 



(56) 
(57) 



IV. COMPARISON BETWEEN THE SCALING SOLUTIONS 



Let us recall the known results concerning the scaling variable s g (|34|). T he beginning of the leading edge expansion 
of the traveling wave solution associated to s g has been calculated in jH, [TB|, and writes, in our notations 



N(L,Y) (x e 



-7cS 9 +0(logY) 



Ai 



V 



with Sn 



L 



1/3 



1/3 



(58) 



(59) 



It is valid for large L and 1 <C s g Y" 1 / 4 . The neglected terms . . . are of order 0(Y 1 / 6 ) when s g — O (Y" 1 / 6 ). 
This section is devoted to the comparison of that solution with the solution (JT]) derived in the previous section. 



A. The saturation scale 



The saturation scale is the typical scale at which occurs the transition between the dilute and the dense regime, i.e. 
between the linear and the nonlinear regime. As the transition is smooth, that definition for the saturation scale is 
ambiguous. For simplicity, we assume that the magnitude of the solution (|56[) is dominated by the exponential factor. 
It leads to define the logarithmic saturation scale L S {Y) by the implicit relation 



Sl (L s (Y),Y)-0(logL) = ( sat , 
with ^ sat some parameter of order one. Then, the usual saturation scale Q S (Y) is defined by 

'Q 2 s(Y)\ 



L S {Y) = log 



A 2 



(60) 



(61) 



Equivalently, it is useful to consider the reciprocal function Y S (L) of L S (Y), i.e. the saturation rapidity at a given 
transverse scale. Using our result ((6]), one gets 



±-*k(DL)i+OQogL)-C 



+ Y . 



(62) 




FIG. 1: Saturation scale. Solid lines: new running coupling scaling solution (162[l ; Dashed lines: running coupling geometric 
scaling solution (|64|l : Dotted lines: fixed coupling geometric scaling solution. Left: Saturation curve in the (L,Y) plane. For 
each scaling solution, the three curves corresponds (from the left to the right) to the values 0.5, 1 and 1.5 for the parameter 
(sat (see equation (|60p ). Right: Derivative of the logarithmic saturation scale &L B /&Y (see equation (JBTJ ) as a function of the 
rapidity Y. 



Note that the Csat dependance, parametrizing the ambiguity in the definition of the saturation scale, is subleading 
compared to the universal terms, even the unknown one O(logL). The subasymptotic Y$ parameter is even less 
relevant. Then we just drop it in the following. Inverting approximately the expansion (|62|) . one finds 

LlOn = ^ + f Di § + * [0(lo S Y) + 2 CsQt ] . (63) 

Finally, the square root of that expression gives 

i 

Ls{Y) = x ( D V^r) 3 + ° {logY) + Csat ' (64) 

which is exactly the expression for the logarithmic saturation scale given by the other running coupling solution 
(|59[) . neglecting the Y Q . Hence, although the two approximate asymptotic solutions of the BK equation with running 
coupling are different, their prediction for the saturation scale agree asymptotically. It is a first evidence that 
those two approximate solutions are approximations of the same exact solution. They agree concerning the first 
subasymptotic correction to the saturation scale, and they would presumably agree for the second one, which is still 
unknown in both cases. Indeed, these first corrections are constrained by the universality property of the asymptotic 
traveling waves solutions. For that reason, the two approximate solutions allow both to extract the leading terms of 
the exact logarithmic saturation scale at running coupling. 

The behavior of the saturation scale determined from the two running coupling scaling variables s g and s n and 
from the corresponding fixed coupling scaling variable are compared in FIG. [I] The saturation line and the derivative 
of L S (Y) (see equation ([HI])) are plotted using the expressions (f6"H) . (f6"2"|) . dropping the 0(logY) or O(logL) terms, 
and taking In = 0. For each solution, the saturation line is plotted for three different values of the parameter Csat- 
As expected pjj], at large Y the saturation line is parabolic for running coupling and linear for fixed coupling. The 
two running coupling solutions agrees concerning the shape, but not concerning the normalization of L S (Y). It comes 
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FIG. 2: N(L, Y) as a function of L for the s g solution (|58|l in dashed line, and for the s n solution (I56p in solid line. From the 
left to the right the different curves corresponds to Y — 2 to Y — 30 by step of 4. The vertical axis is in logarithmic (base 10) 
scale. The two plots just differ by the range they cover. 



partly from the dropped terms, and partly from the ambiguity in the definition of the saturation scale. The parameter 
(sat affects indeed more the the normalization of L S (Y) than its shape. The derivative of L S (Y) corresponds in the 
saturation models (e.g. [33| and many following models) to the parameter usually called A. One remarks that the 
fixed coupling solution is not consistent with the fitted value A ~ 0.3, A being stable in the rapidity range covered by 
HERA. By contrast, such low and stable value can be provided in the running coupling case, as noticed in [T3|, but 
seems to require a large rapidity shift Yq, which may be considered as a signal of important NLL corrections. The 
shape of the saturation scale is indeed less sensitive than the additionnal Yq parameter to higher orders. 



To compare the solutions ([55]) and ([55]) further, let us now consider the shape of N(L,Y) obtained in both cases. 
They are plotted on FIG.0 down to the vicinity of the saturation scale, but not beyond, as we have explicit solutions 
only in the linear regime. The values chosen for the parameters are A = 1, B = 4 and Yq — 0. The left hand 
plot corresponds to a range more relevant for phenomenology, but the right hand plot is interesting concerning the 
mathematical meaning of these two approximate solutions. They are both derived with the traveling wave method, as 
the beginning of an expansion in the leading edge region. As the two solutions coincide over a part of the linear region, 
it is obvious that they correspond to different expansions of the same exact solution of i.e. the exact pulled front 
solution selected dynamically. The range over which the two solutions agree is growing with Y and drifting towards 
the tail of N, as expected for the leading edge region. Hence, the two approximate solutions {55J) and (|56p are in 
agreement concerning all the features constrained by the universality property of the pulled front solutions, i.e. the 
evolution of the saturation scale and the shape of the front in the leading edge region. They give presumably an 
accurate description of these features, at least in the large L limit (at s g or s n fixed). However they differ concerning 
other features, such as the shape of the front around the saturation scale (i.e. in the front interior) or in the tail 
forward the leading edge. For exemple, asymptotic freedom effects are stronger for the s n solution (|56I57[) than for 
the s g solution (158|59|) . as the tail of the front evolves slower than the bulk in the former one. The left hand plot 
shows that near the saturation scale, the front given by the s n solution is significantly steeper than the one given by 
the s g solution. One can study the shape of the front more precisely by calculating an effective anomalous dimension 



B. The shape of the front 




FIG. 3: Effective anomalous dimension 7 e //(L, Y) (as defined by equation ([65])) as a function of L, calculated for the solutions 
s g solution (left panel) and s n solution (right panel). From the left to the right the different curves corresponds in both cases 
to Y — 2 to Y — 10 by step of 2. The horizontal line stands for the value of the critical saturation exponent 7 C ~ 0.63. 



defined as 

%ff (L,Y) = -d L \ogN(L,Y). (65) 

The behavior of j e ff(L, Y) for the two running coupling approximate solutions is compared on FIG. [31 They give 
very different results. Moreover, both differ from the naive expectation j e ff(L,Y) ~ 7 C . The factor containing 
the Airy fonction indeed contributes to the shape of the front in the leading edge, and the nontrivial structure of 
the scaling variable s n gives also a contribution to j e ff(^L,Y) even in the front interior. Therefore, the theoretical 
solutions (|58p and (|56p (and the fixed coupling one) involve the exponent j c , but it does not mean that the anomalous 
dimension of these theoretical solutions is indeed 7 C , especially at finite L and Y. For completeness, it is useful to 
recall the behavior of the universal travelling wave solution of the fixed coupling BK equation ([3J). This solution 
and the corresponding 7 e //(L, Y) are plotted on FIG. [U Its front has a shape very different from the one of the s n 
solution, but rather similar to the one of the s g solution. 

These observations allow to address a puzzle raised by the numerical simulations of the BK equation with running 
coupling [1^, [3(| • It was noticed that numerical simulations leads to a wave front much steeper at running coupling 
than at fixed coupling, whereas one expect a similar wave front from the theoretical fixed coupling solution and s g 
running coupling solution (|58| . We have shown that the latter is only an approximate solution of (j4|), and that ([56]) 
is another one, within a complementary approximation. Hence the exact pulled front solution of ^ lies probably 
between these two approximate solutions, and can then have a stronger "f e ff(L,Y) than the fixed coupling solution 
in the front interior. 



V. PHENOMENOLOGY 

In this section, we discuss the phenomenological validity of the new scaling, using the results of Ref.[32|. The 



quality factor method 31| allows to test and quantify scaling behaviors of experimental data. For a given data set, 
the higher is the quality factor value, the better is the scaling behavior. This method is used in [32] to compare scaling 
behaviors at the cross section level. Indeed, geometric scaling at the cross section level comes from geometric scaling 
properties of the dipole-target amplitude through dipole factorization and the replacement k 2 1— > Q 2 in the scaling 
variable. Note that such property is only approximate for the new scaling variable s n . 
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FIG. 4: Left: N(L, Y) as a function of L for the fixed coupling solution fltj [Isl] in dashed line, and for the s n solution (|56[1 in 
solid line. From the left to the right the different curves corresponds to Y = 2 to Y = 30 (or Y — 26 for fixed coupling) by 
step of 4. The vertical axis is in logarithmic (base 10) scale. Right: Effective anomalous dimension y e ff(L,Y) (as defined by 
equation (|65p) as a function of L, calculated for the fixed solution. From the left to the right the different curves corresponds 
in both cases to Y = 2 to Y = 10 by step of 2. The horizontal line stands for the value of the critical saturation exponent 
7 C ~ 0.63. 



Starting from the theoretical scaling variables st (|2J), s g (|3~4"|) and s n ([55]) , we shall use the following versions at the 
cross section level 

A Y , (66) 



A y/Y-Yo, (67) 

A is taken as a free parameter in each case. The rapidity shift Yq is either fixed to zero or treated as a free parameter. 
Note that such shift cannot affect the quality of the scaling with the 77 variable. The value of A has an influence on 
the scaling quality only in the r n scaling case. A is either fixed to Aqcd — 0.2 GeV, or treated as a free parameter. 

Using the quality factor method on the inclusive D1S cross section data from HI, ZEUS, E665 and NMC experiments 
(see [32| for a complete discussion of data selection and references) in the range 3 GcV< Q 2 < 150 GeV and x < 10 -2 , 
we have found the following results [32j. The optimal quality factor from one scaling variable to another, when A is 
the only free parameter, do not differ more than 4%. The data thus scale equally well with the three variables. The 
results for the new variable are displayed in F1GJ5J On the right-hand plot, the data shows a nice scaling property 
with T n , which is effectively as good as with r/ [IlTl3lT| or as with r ff [3lT |. On the left-hand plot, one remarks a sharp 
peak in the quality factor, showing that the parameter A is determined unambiguously. That peak is remarquably 
stable when the study is done on some subsets of the data, like Hl+ZEUS data only, or HI only, or ZEUS only. It 
illustrates the robustness of the quality factor fit in that case. 

When introducing the other parameters Yq and A in the quality factor fit, we find that for the running coupling 
geometric scaling variable r g , the optimal value for Yq remains compatible with zero. By contrast, for the new scaling 
with t„, letting A and Yq free improves significantly the quality factor from QF — 1.69 to QF — 1.82 (in 10~ 3 units). 
The optimal values of the additional parameters are Yq = —1.2 and A = 0.3 GeV, which are reasonably moderate 
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FIG. 5: Phenomenological results for the new running coupling scaling variable r n (|68[) . with Yo = and Aqcd = 0.2 GeV. 
Left: Normalized quality factor as a function of the parameter A for different data sets. Right: Inclusive DIS cross section as a 
function of the scaling variable r n (Q 2 ,Y). 



shits. Hence, the prefered scaling behavior found in the inclusive DIS cross section data is obtained using the new 
scaling variable t„. This new scaling works as well for DVCS, vector meson production, and diffractive data [32| . 

VI. CONCLUSION AND OUTLOOK 

We have shown that the Balitsky-Kovchegov equation with running coupling has no exact scaling solution, contrary 
to the fixed coupling case. However, that equation is compatible with two types of approximate scaling solutions, 
which we have analyzed. The first type of scaling behavior boils down to the running coupling geometric scaling ((5]), 
already known while the second one §§§ is new, and goes beyond the notion of geometric scaling. 

Using the traveling wave method, we have derived the asymptotic traveling wave solution (|7I8[) associated with 
the new scaling variable ©. Both running coupling traveling wave solutions are shown to be compatible concerning 
the evolution of the saturation scale, and the shape of the leading edge of the pulled front, while they differ in the 
front interior. It means that the asymptotic traveling wave solutions associated with the two scaling variables are 
both approximations of a same exact solution. The asymptotic shape of the leading edge and the evolution of the 
saturation scale are universal features of the exact solution, and the traveling wave method allows to reach them, 
within both approximations, because both approximate scalings become exact asymptotically. 

By contrast, other features such as the shape of the front closer to the saturation scale arc more difficult to 
calculate within the traveling wave framework. One finds that the front is steeper (see FIGH|) for the new solution 
([7|) than for the geometric scaling solution (|58|) . As the two approximate solutions rely on similar but distinct 
approximations, one could estimate their uncertainty by their difference. The exact solution is presumably lying 
between the two approximate solutions. Hence, it should have a steeper front than the running coupling geometric 
scaling solution, and than the fixed coupling solution, as found in numerical simulations |29l. 1301] - 

In the fixed coupling (and mean field) case, the main geometric scaling violations are the universal ones resulting 
from the convergence towards the asymptotic behavior. In the running coupling case, the dipole target amplitude 
shows simultaneously the two scaling behaviors (0 and ©. As these two scalings are compatible only in a large L 



and large Y limit, that property should bring violations of both scalings at nonasymptotic energies, in addition the 
universal scaling violations. Those stronger scaling violations at running coupling, and the fact that the new scaling is 
slightly favoured by the DIS data (3^ |. explain why the running coupling geometric scaling is not significantly better 
than the fixed coupling geometric scaling at HERA energies. 

The NLL and higher order terms are quantitatively important at HERA energies [l?} • The phenomenological 
determination of the scaling parameters done in [32j and in section [V] may be useful to understand the nonasymptotic 
effects of those contributions. 

Let us briefly discuss our various theoretical approximations. Following j2(jj |, we have neglected Pomeron loops 
effects, which seem irrelevant in dense-dilute collisions in the running coupling case. Our result (|7|8p is free from 
NLL corrections. They appear indeed at the next order in the leading edge expansion |7| only. However, NLL and 
higher orders are phcnomenologically important, as mentionned previously. Our choice for the running coup ling scale 
differs from the best ones [HI, [3] , but this difference is only marginally relevant in numerical simulations [30] . The 
dipole or parton correlations present in the Balitsky-JIMWLK equations, and the off-diagonal momentum transfcrt 
dependance, do not appear in the linear part of the equation. Hence, the universality of the pulled front solutions 
guaranties that these correlations cannot modify our result (|7I8|) . Note that the effect of all those approximations is 
obviously the same for the previous approximate solution (|58I59|) . 

As an outlook, it would be interesting, on the theoretical side, to explore the tail of the front. The two scaling 
behaviors are no more relevant when s n © or s g ([5]) are of order L, i.e. in the forward region of the front, where the 
initial condition still propagates. Hence, the crossover between the pulled front and the forward region seems more 
complex than in the fixed coupling case. It should be also important to calculate the next order of the leading edge 
expansion of the two running coupling solutions, including NLL effects. Moreover, the variable s n may be useful to 
study the solutions of the (LL or NLL) BFKL equation with running coupling. 

On the phenomenological side, a systematic study [321 ] of scaling properties of inclusive and exclusive DIS cross 
sections has been conducted simultaneously to the present paper. A further step would be to use our results to 
build a running coupling saturation model, in a similar way as the Iancu-Itakura-Munier model [3~ij for the fixed 
coupling case. Most available phenomenological models of saturation are built at fixed coupling. Although successful 
to describe HERA data, they may fail to give reasonable predictions for the LHC, due to the different evolution of 
the saturation scale. This justifies the need for a running coupling saturation model, e.g. based on the new scaling 
variable. 
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APPENDIX A: DERIVATION OF THE LEADING-EDGE OF THE NEW ASYMPTOTIC SOLUTION 

We finish in this appendix the calculation of section HTT1 starting from Eq. 



= x'ilc) 



fLE(Sn,L)+J2^X {P) (lc) 



Or 



IcSr, 



As discussed in section |HI| one can factorize f^E as 

f LE {s n ,L)~e-^ Ll,Z G(z) 

z being 

3/3, 



f LE (s„,L). (Al) 



(A2) 



Zy3 L3 



(A3) 



with the shifted scaling variable 



Inserting (|A2[) in the equation (| Al[) written in z and L variables, i.e. doing the replacements 



d s 



one gets 
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G(z). (A4) 



G(z) is determined by the terms of order L 2 / 3 , for z = 0(1). Hence we have 

X '( 7c ) 7c (z - 4/3) G(*) = Jx"(7c) G"(z). 



(A5) 



The following terms in (|A4[) are of order L^ 1 . The subleading term in Jle is thus suppressed by a factor L~ 1//3 , 
compared to G(z). A factor L s appears when calculating the subleading term in /lb, in the same way as the factor 
exp(— 37C/3-L 1 / 3 ) appears when calculating the leading term. Our choice is to single out an L 1 / 3 factor from this L s , 
and to write the remaining part as a 0(logL) correction to the scaling variable. Finally, Jle writes 



f LE ( S , L) = e-WLV'+OQegL) L l/3 J G(z) + L -k H{z ) + O 



(A6) 



for L large, and z — 0(1). The function H(z) would be determined at the next order, as the 0(logL) correction to 
the scaling variable. 

Let us now calculate G{z) and j3 from the equation (|A5|) . Up to an affine change of the variable z, this equation 
can be mapped into the Airy equation. Hence G{z) writes 



G(z) = Gi Ai 



4/3 



Di 



G 2 Bi 
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Di 



(A7) 



where Ai and Bi are the Airy functions, C\ and G2 are two integration constants and D, playing the role of a diffusion 
coefficient, is defined by 
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X"(7c) X"(7c) 



2 7cX'(7c) 2 X (7c) 



(A8) 



Our calculations are valid only if s n <C y/L. A priori, it allows s n to be large compared to (-DL) 1 / 3 , and then z to 
be arbitrary large for L — > 00. Thus, G(z) has to be regular for z — > 00, which gives G2 = 0. The other boundary of 
the interval of validity of the solution (|A6jl is located at z ~ 0, and corresponds to the transition to the saturation 
regime. For 1 <C s n L 1 / 3 , the solutions (|A6|1 and (|47|) are both valid, and then we have to match them. We 
have determined that the scaling is modified by the leading edge dynamics from s n to s n . And in the front interior, 
the difference between s n and s n is a subleading effect for the solution (|4T]) . Hence, we can match the leading edge 
solution 



N(s n ,L) 



_ - lc s n +0{\ogL) 
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with the front interior solution with s n replaced by s n 

N (s n ) = A (s„ + B) 
For 1 <C s„ -C L 1 / 3 , (|A9|1 can be expanded as 



N{ Sn ,L) = e -7cS„+0(logL) L i 
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(All) 



Hence, one finds that — 4/3/Z) 1 / 3 has to be a zero of the Airy function Ai. Requiring that N(s n ,L) is positive for 
z > selects the rightmost zero £1 ~ —2.338. Hence, one has 



The matching of the terms proportional to s n gives 



£>l/3 
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The constant terms give then H(0) = AB. However, one can replace s n by s n + B in the argument of the Airy 
function Ai in (|A9|) . It corresponds to a choice of the repartition of the higher order terms, very similar to the choice 
of the L 1 / 3 in front of the solution (|A9j) . Using this prescription, one has H(0) = 0, and then the contribution of the 
higher orders is reduced, at least in the neighborhood of z = 0. 

Finally, the universal asymptotic traveling wave solution of the BK equation with running coupling associated with 
the scaling (j3"5)) writes, for large L and 1 -C s n <C yL 



N(L,Y) 
with s n 
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Note that we have not determined the parameters A and B. However, they are not completely free, and they should 
be constrained by matching the expression (|47p with numerical, or analytical if possible, solution of the nonlinear 
equation (f3"8"]) with v — v c . Note also that the prefactor (DL) 1 / 3 is arbitrary, as it is part of the still unknown O(logL) 
term in the exponent. We have written it just for convenience, in order to cancel in the large L limit the (DL) 1 / 3 
denominator present in the Airy function. The same is true for the Y" 1 / 6 , or i 1 / 3 , prefactor in the solution associated 
the variable s„ in [la 11811. 
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